Theory of spin Hall effect 
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An extension of Drude model is proposed that accounts for spin and spin-orbit interaction of 
charge carriers. Spin currents appear due to combined action of the external electric field, crystal 
field and scattering of charge carriers. The expression for spin Hall conductivity is derived for metals 
and semiconductors that is independent of the scattering mechanism. In cubic metals, spin Hall 
conductivity o~ s and charge conductivity a c are related through a s = [2-Kh/(3rnc 2 )]a 2 with m being 
the bare electron mass. Theoretically computed value is in agreement with experiment. 

PACS numbers: 72.25.-b, 72.10.-d, 71.70.Ej 
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It has been a common knowledge in atomic physics 
that due to spin-orbit interaction the spatial separation of 
electrons with different spin projections can be achieved 
through scattering of an unpolarized electron beam by 
an unpolarized target Dyakonov and Perel were the 
first to notice that in the presence of the electric cur- 
rent the scattering of charge carriers by impurities in a 
semiconductor must lead to a similar effect It was 
subsequently called the spin Hall effect 0] and observed 
in semiconductors [J, [5| and metals [6| . A number of mi- 
croscopic models have been developed that explain spa- 
tial separation of spin polarizations by various "extrinsic" 
(due to impurities) and "intrinsic" (impurity-free) mech- 
anisms, see for review Ref.0 While these models provide 
valuable insight into microscopic origin of the spin Hall 
effect, they are lacking universality of, e.g., Drude model 
of charge conductivity [§]. The Drude model, in spite of 
being classical in nature, has been very powerful in de- 
scribing dc and ac conductivity and its temperature de- 
pendence. It also gives the accurate value of the Hall co- 
efficient by catching correctly the orbital motion of charge 
carriers in the presence of the magnetic field. The power 
of the Drude model resides in the fact that it expresses 
conductivity, an = e 2 n,T/m, via charge e, concentration 
n, mass m, and relaxation time r of charge carriers re- 
gardless of the scattering mechanism. Same parameters 
enter expressions describing experiments other than the 
Ohm's law, e.g., n and the sign of e can be extracted from 
measurements of the Hall coefficient Rh = — (nec) -1 , r 
can be extracted from measurements of the frequency de- 
pendence of the impedance, and m can be extracted from 
measurements of the cyclotron resonance. This allows 
one to test theoretical concepts of charge conductivity 
regardless of the degree of accuracy with which one can 
compute parameters entering a d- 

In this Letter we will try to develop a similar ap- 
proach to the spin Hall conductivity. We will take the 
Drude model a little further by incorporating spin and 
spin-orbit interaction into the dynamics of charge carri- 
ers. We will argue that such a straightforward extension 



of the Drude model allows one to obtain universal ex- 
pression for spin Hall conductivity that is independent 
of the scattering mechanism. The spin Hall effect ap- 
pears due to combined action of the external electric field, 
quadrupole crystal electric field and scattering of charge 
carriers. Same as for charge conductivity, all details of 
the scattering mechanism are absorbed into the momen- 
tum relaxation time r. We will show that this crude 
model provides correct values of spin Hall conductivity 
in both metals and semiconductors. 

For certainty we will speak about electrons but the 
model will equally apply to holes. The non-relativistic 
limit of the Dirac Hamiltonian for spin-1/2 particle is Q 
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Here U(r) represents the action of microscopic electric 
field on charge carriers. It incorporates effects of electro- 
static crystal potential $ , potential due to imperfections 
of the crystal lattice <&i, and external potential $ e , 



U(r) = e$ (r) + e$. ( (r) + e$ e (r) . 



(2) 



The last term in Eq. |T]) is the spin-orbit interaction, with 
<r being Pauli matrices. 

Hamiltonian mechanics for canonically conjugated 
variables p and r is described by the equations 
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(These equations can also be derived from quantum- 
mechanical relations: ihr = [r, H], ihp — [p,H], 
[fitPj] — ifiSij.) From Eq. (J3J) one has 
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The second term in the expression for the momentum, 
Eq. ([5]), that is proportional to the cross product of the 
electron magnetic moment and the electric field, is the so- 
called "hidden mechanical momentum" associated with 
the momentum of the electromagnetic field, see e.g. Ref. 
Hol and references therein. 

Substitution of Eqs. © and © into Eq. fU) yields 
the following form of the second Newton's law for charge 
carriers: 
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dr 



F CT (r,r). 



(7) 



where the spin-dependent force is given by 
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Here we neglected the term proportional to 1/c 4 that 
exceeds the accuracy of Eq. (fTJ) and took into account 
that r in the second line of Eq. ([5]) , that originates from 
the expression for p in Eq. ([5]), should not be differen- 
tiated on r because p and r in Eqs. © and ^ are in- 
dependent canonically conjugated variables. Note that 
the force in Eq. ([5]) is equivalent to the Lorentz force, 
F CT = (e/c)(r x B ff ), acting on a particle of charge e in 
the magnetic field 
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(a x E tot ) , (9) 



with E t ot being the total electric field, eE tot = ~dll/dr. 
One can trace this force to the fact that with an accuracy 
to c~ 2 the Hamiltonian fTJ) can be written as [ll[ 
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U(r). 
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Crystal field creates a non-zero average of this force any 
time the charge carriers have a non-zero drift velocity, 
(r) ^ 0. It is this effective Lorentz force that is re- 
sponsible for the spin Hall effect. This our conclusion 
is similar to the conclusion of Hirsch [l2j who studied 
the force exerted on a line of moving magnetic dipolcs 
by the electrostatic field of charges arranged in a cubic 
lattice. For a cubic lattice (see below) our result for the 
effective Lorentz force coincides up to a factor of 2 with 
the result obtained by Hirsch. As in our approach, Hirsch 
found that the effective Lorentz force is generated by the 
second derivative of the crystal field. In computing this 
force he replaced moving magnetic dipoles with station- 
ary electric dipoles that produce an equivalent electric- 
field. This eliminated the hidden momentum responsi- 
ble for the first term in the second line of Eq. JSJ • The 



absence of this term in Hirsch's model accounts for the 
above mentioned difference by a factor of 2. 

In the spirit of the Drude model we shall now add to 
Eq. ([7]) the drag force — mi jr. Then Eq. ([7]) becomes 
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Here we assume that to the first approximation the ve- 
locity relaxation of charge carriers is independent of their 
spin, that is, r is independent of <x. Since the relaxation 
is due to imperfections of the crystal lattice, in order not 
to count their effect twice, we should now think of U (r) in 
Eqs. (fTT|) and ([8|) as a sum of the ideal periodic potential 
of the crystal lattice, e$0j and the potential produced by 
the externally applied voltage, e<E> e . Due to relativistic 
smallness of the spin-dependent force ([5]) one can treat 
F(j(r, r) in Eq. (flT]) as a perturbation. Then the solu- 
tion of Eq. (jTTj) can be written in the form r = ro + rj., 
where v\ is a small spin-dependent part of the velocity 
proportional to c~ 2 . In the presence of a constant exter- 
nal electric field E = —d$ e /dr, with the linear accuracy 
on E, one obtains from Eqs. (fTTj) and (JSJ) 
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In deriving Eq. (|13|) we have made an assumption that 

<F CT (r ,r )) = - <r X B ff (r )} = - (r ) x (B ff > . (14) 
c c 

Some justification of this assumption is provided by the 
following argument. By symmetry, the only reason for 
(F CT ) to be different from zero would be (r ) ^ 0. Con- 
sequently, (F CT ) should be first order on (r ). Being per- 
pendicular to the velocity, the force F a — (e/c)(r x B CT ) 
does not do mechanical work on the charge. Neither 
should (F CT ) with respect to (r ), rendering the form 
(F CT ) = (e/c)((r ) x & e ff)- Since the trajectory of 
the particle ro(t) does not have strong correlation with 
the quadrupole component of the crystal electric field, 
the above factorization of the average with the choice 
B e // = (Bo-) should not deviate strongly from the exact 
average. 

The right hand side of Eq. (flU)) contains the volume 
average of V;Vj<I>o( r )- In what follows we will study 
the spin Hall effect in a cubic lattice. Generalization to 
other lattices is straightforward and will be considered 
elsewhere. The case of a cubic lattice, besides simplic- 
ity, is interesting because experiments performed to date 
have been done in cubic semiconductors and in aluminum 
that is also cubic 0, [H, Q • In the cubic case the only in- 
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variant permitted by symmetry is 

d 2 $ \ _ M 
dridrj I 13 ' 
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where A is a constant to be determined later. With the 
help of Eq. (|15[) one obtains from Eq. (Qj: 
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Let us now introduce the vector of spin polarization 
of the electron fluid, £ = (a), which absolute value lies 
between and 1, 



n + — 7i- 



(17) 



Here n± are concentrations of charge carriers with spins 
parallel and antiparallel to £ respectively, with n++n_ = 
n being the total concentration of charges carrying the 
electric current. For a ferromagnet, £ ^ is the equi- 
librium state of charge carriers below the Curie temper- 
ature, while for a small non-magnetic conductor a sig- 
nificant value of £ can be achieved through injection of 
spin-polarized charge carriers from a magnetic conduc- 
tor, see for review Ref. 13 The density matrix of the 
charge carriers in the spin space (normalized to their to- 
tal concentration n) can be written as 



TV 



The electric current is 



(18) 



j = e(Arr) = e(iV(ro + r 1 )). (19) 

Substituting here equations (|12|) and (fT6|) one obtains 

j = (7 c E + £r,KxE], (20) 

where charge conductivity and spin Hall conductivity are 
given by 



e 2 nT 



2m 3 c 2 

An interesting observation is that the ratio 
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is independent of the concentration of charge carriers. 
This result also follows from microscopic models of spin 
Hall effect 0, EH. It explains why this ratio has the 
same order of magnitude in metals and semiconductors 
0, [H] despite very different concentrations of charge car- 
riers. Note that according to Eq. (|2"3")l the temperature 



dependence of the ratio of spin Hall and charge conduc- 
tivities is determined by the temperature dependence of 
the relaxation time r. Consequently, in metals one should 
observe monotonic temperature dependence of spin Hall 
conductivity, cr s (T) oc a 2 {T), while in semiconductors 
a s (T) oc n(T)r 2 (T) may exhibit maximum on tempera- 
ture [161 ]. 

Eq. (|13j) shows that quantitative analysis of the spin 
Hall effect requires knowledge of the volume average of 
the quadrupole component of the crystal electric field, 
(d 2 $o/ dridrj). For a specific crystal it can be computed 
by methods of density functional theory. For a cubic 
crystal (d 2 $o/dridrj) reduces to a single constant A, see 
Eq. (fT5|) . Here we will provide a simple estimate of this 
constant. In a metal the electrostatic potential created 
by the crystal lattice of positively charged ions satisfies 
Maxwell equation 



V 2 $ = -4np(x, y, z) 



(24) 



where p(x,y,z) is the local charge density of ions. The 
volume average of this equation is 
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dy 2 
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dz 2 



= 3A = -4Tr(p(x,y,z)) 
(25) 

Due to the periodicity of the crystal field, it is sufficient 
to compute this average over the unit cell of the crystal. 
This gives 
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where — Ze and hq are the charge and the concentration 
of ions respectively. The sign in (p) = —ZeriQ is deter- 
mined by our choice of the negative charge e of the charge 
carriers. In the accurate quantum-mechanical calcula- 
tion of the parameter A produced by the distribution of 
charges in a cubic lattice, Z should come out as a number 
of order unity when hq is chosen as the inverse volume 
of the unit cell. This should be true for both metals and 
semiconductors . 

Substitution of Eq. ([H into Eq. J23D gives 



(27) 
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At n - 10 22 cm- 3 (10 23 c m - 3 ) and r - l(T 13 s (10~ 14 s) 
this ratio is of order 10 which is in agreement with re- 
cent experimental findings in metals and semiconductors 
0, [||]. According to our model such a value of cr s /a c is 
not universal though. In pure conductors, when scatter- 
ing of charge carriers is dominated by phonons, large r 
at low temperature can provide much greater values of 
the ratio of spin Hall and charge conductivities. 

For the absolute value of spin Hall conductivity one 
obtains from Eq. ([22| 
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In metals, when Z is identified with the ion valence, one 
can replace Zuq with n and write Eq. (|28[) as 



2nh 
3mc 2 



(29) 



This relation between spin Hall and charge conductivities 
permits comparison with experimental data without any 
fitting parameters. Such a comparison with the data on 
aluminum [6j] is shown in Tabled] (Note that conductiv- 



Conductivity 


Experiment (Q 1 m x ) 


Theory (fi^m -1 ) 


a c (12nm) 


1.05 x 10' 




Us (12nm) 


(3.4 ±0.6) x 10 J 


2.6 x 10 d 


a c (25nm) 


1.7 x 10' 




g s (25nm) 


(2.7 ±0.6) x 10* 


6.9 x 10 a 



TABLE I: Comparison of Eq. (I29p with experimental data }6i] 
on charge and spin Hall conductivities in Al strips of 12-nm 
and 25-nm thickness. 

ity in Gauss units that have been used to derive Eq. (f2"9")l . 
equals conductivity in units of f^m -1 times 9 x 10 9 .) 
Given the crudeness of the model the agreement by or- 
der of magnitude between theory and experiment is quite 
remarkable. Note that in the spirit of Drude model one 
can also make a straightforward extension of our model 
to the ac spin Hall conductivity. 

Throughout our derivations we did not distinguished 
between the effective mass and the bare mass of charge 
carriers. To account for their difference we notice that 
the spin-orbit term in Eq. ([T]) can be interpreted as the 
Zeeman interaction of the electron spin with the effec- 
tive magnetic field that is generated by the electrostatic 
field in the moving frame of the electron. Consequently, 
one of the masses in S/(4m 2 c 2 ) in front of the spin-orbit 
term must be the bare electron mass that enters the Bohr 
magneton, while the other m is due to orbital motion. In 
a solid, one should replace m 2 in the spin-orbit term with 
mm* , where m* is the effective electron mass. Also in 
the first term of Eq. (Q]) m should be replaced by m*. 
After these replacements are made and followed down to 



Eq. 



gets absorbed into a c while 1/m in front 



of a c can be traced to the expression for the Bohr mag- 
neton, \xb = \e\h/(2mc). Thus the coefficient in front of 
a 2 in Eq. (f2"9")l contains the bare electron mass. We used 
this fact for computing numerical values of a s presented 
in Table [J 

The theory of spin Hall effect developed in this Letter 
should also apply to the anomalous Hall effect observed 
in ferromagnets. According to Eq. (|16[) the weak flow of 
charge carriers with opposite spin polarizations in the op- 
posite (perpendicular to E) directions occurs regardless 
of the average polarization £ given by Eq. (jTTJ). When 
the time of spin relaxation of charge carriers is sufficiently 



long, this will lead to the spin polarization of the bound- 
aries of a small conductor even when £ = 0. On the 
contrary, the spin Hall current given by the second term 
in Eq. (|20|) is proportional to £. In non-magnetic metals 
and semiconductors, large spin polarization is not easy 
to achieve but in a ferromagnetic metal £ can be of order 
unity, resulting in the noticeable current, cr s [£ x E], that 
is normal to E and normal to the magnetization. 

In Conclusion, we have presented a simple extension of 
Drude model of conductivity to the case of charge carri- 
ers that have spin and spin-orbit interaction. Spin Hall 
effect appears naturally in such a theory as a combined 
action of the external voltage, crystal electric field, and 
scattering of charge carriers. The expression for spin Hall 
conductivity is independent of the mechanism of scatter- 
ing. Theoretical values of spin Hall conductivity com- 
puted within such a model are in agreement with exper- 
imental data on metals and semiconductors. The model 
can, therefore, serve as a simple physical picture of spin 
Hall effect. 
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